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Motivated by neutron scattering experiments in tlie liigh-Jc cuprates, we propose tlie two-component spin- 
fermion model as a minimal phenomenological model which has both local spins and itinerant fermions as 
independent degrees of freedom. Our calculations of the dynamic spin correlation function provide a successful 
description of the puzzling neutron experiment data and show that: (1) the upward dispersion branch of magnetic 
excitations is mostly due to the local spin excitations; (2) the downward dispersion branch is from collective 
particle-hole excitations of fermions; and (3) the resonance mode is a mixture of both degrees of freedom. 
Using the same model with the same set of parameters we calculated the renormalized quasiparticle dispersion 
and successfully reproduced one of the key features of the angle resolved photoemission spectroscopy (ARPES) 
experiments, i.e., the high energy kink structure in the fermion quasiparticle dispersion, hence further support 
^ , the two component spin-fermion phenomenology. 
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I. INTRODUCTION 



The study of spin dynamics has been a key research interest since the discovery of the high-T^ superconductor because it is 
expected that the spin correlation holds crucial information for the mechanism of the high-71 superconductivity (HTS). For long 
the two main observations in the neutron scattering experiments of high-F^ cuprates (HTC) are (1) the incommensurate (IC) peaks 
"j^ at low energy or at quasielastic excitations^'^ and (2) the so-called resonance peak at commensurate wave vector at relatively high 
energy (30 ^ 50 meV)^"^. In early experiments the IC peaks were observed only in the deeply underdoped lanthanum cuprate 
I and the resonance mode was the hallmark of the fully doped two layer yttrium cuprate. However, later experiments reveal 
that both features appear in both groups of cuprate compounds, although sensitively depending on doping level. More recently 
C with inelastic neutron scattering (INS) experiemnts^^' with high precision, an unifying form of the magnetic excitations in the 
cuprate superconductors has emerged as "hourglass" shape of excitations around the wave vector (1/2,1/2) (hereafter in units of 
2%jd), in which the low energy IC excitations form the downward dispersion branch and the high energy IC excitations form the 
upward dispersion branch, and the two branches of excitations merge at the commensurate momentum Q=(l/2,l/2) and at the 
resonance frequency Q-res- 

It is a pressing question to understand the origin of this "hourglass" shape excitations. Theoretical proposals up to now can 
be classified into two groups: (1) theories based on the spin dynamics in the presence of stripes'^"^^, and (2) Fermi liquid type 
theories of itinerant fermions The key idea of the first group of theories is that the stripes formed by doping in the two 
dimensional Cu-0 plane splits the commensurate spin wave excitations into two IC branches at the wave vectors (1/2 ±5, 1/2) 
or at their symmetry rotated positions by x < — > y depending on the directions of the stripes. The dispersions from each branch 
of the two IC modulation cross at the commensurate wave vector (1/2,1/2) at a higher energy, which is then identified as the 
resonance mode. This picture provides a qualitative explanation to the hourglass dispersion and the resonance mode. However, 
this type of theories has difficulty to be extended to the higher doping regime where the presence and the nature of the stripes 
is questionable. The second group of proposals are itinerant fermion theories with interaction^^— . In this type of theories, 
the resonance mode and the downward dispersion can be obtained, but the upward dispersion branch is not yet satisfactorily 
reproduced. 

I In this paper, we propose a two component spin-fermion model^''^^ as a minimal phenomenological model to provide a 
natural and unifying explanation of the above mentioned neutron experiments of HTC. In this phenomenological model, the 
minimal set of low energy degrees of freedom are the spin wave excitations of local spins and the continuum particle-hole 
excitations of fermions. A similar phenomenological theory is also known as one component spin-fermion model and has been 
intensively studied by Pines and coworkers^^. The main difference of the two component model from the one component one 
is the introduction of the spin wave excitations directly from the local spins in addition to the usual collective spin density 
excitations from fermions. In this paper, we show that the presence of the local spin fluctuations is essentially proven by the 
INS and angle resolved photoemission spectroscopy (ARPES) experiments of HTC, therefore supporting the two component 
spin-fermion model as a minimal phenomenological model of HTS. 
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II. FORMALISM 

In a mixed momentum and real-space representation the two component spin-fermion phenomenology Hamiltonian is written 

as 

// = ^4(k)e(k)ca(k)+ ^ ^S(r)-4(r)a„pcp(r)+//s(S(r)), (1) 

k,a r.a.p 

where the first term is the fermionic kinetic energy and the second term describes the coupling between local spins S(r) and 
the spin density of the conduction electrons s(r) ~ c^(r)aapC|3(r). The last term Hs{S{r)) represents an effective low-energy 
Hamiltonian for the local spins S(r). Instead of specifying Hs and solving it, we assumed a phenomenological Ansatz of the 
bare (before coupling to the fermions) local spin dynamics < S^S-^^ >n= Xo.s{(l,^) with a short range AFM correlation, which 
has the general form as follows^^. 

Xa^(q,i2) = Xas(Q^O) • [1 + ^'|q - Q|' - ^V^g] , (2) 

where Q the 2D AFM ordering vector, and the spin gap energy Asg and the magnetic correlation length t, combine to give the 
spin wave velocity — Asg ■ ^ which can be determined by direct measurement^^. 

The key difference of the Hamiltonian Eq.(l) from the one component Hamiltoniar^ is the definition and meaning of the spin 
fields S(r). We assumed that S(r) is the local spin degrees of freedom (d.o.f.) besides and independent from fermions Ca(r) and 
related itinerant spin density s(r) = Ca(r)Oor|3Cp(r). On the other hand, in one component spin-fermion model^^, the collective 
spin fields is defined as the itinerant spin density operator made of fermions, hence fundamentally linked to the fermions, and 
there is no concept of the local spins. Another important distinct feature of our two component phenomenology is the local spin 
dynamics defined in Eq. This form of the local spin correlation function with a short range AFM order should be valid not 
only near Q but also for the entire BZ of q - if the precise form of the spin-wave dispersion is ignored - because at high energies 
the local spin dynamics becomes less sensitive to the long range order or short range order. However, the similar form of the 
itinerant spin correlation function assumed in the one component model^^ is valid only in a narrow region of q around Q by 
definition. 

Microscopic justification of the above two component model, starting, for example, from Hubbard or t-J model, is the heart of 
problem of HTC for the last twenty years or so. We can only sketch here the underlying idea for our phenomenology. Starting 
from a Hubbard model, for example, dynamic mean field theory (DMFT)^^ demonstrated that the key consequence of the strong 
correlation of large U Coulomb interaction is to split the electron spectral density into two parts: one near the Fermi level - the 
itinerant one, and the other at the lower and upper Hubbard bands far below and above the Fermi level - hence the localized one. 
Therefore this splitting of one bare electron spectral density into the itinerant and localized parts is not a new observation but 
has already had a solid theoretical justification. 

A new step in our phenomenology is to propose that the localized spectral density far away from the Fermi level is not dormant 
for the low energy physics. In the framework of the DMFT, once the coherent band is formed at the Fermi level in addition to the 
upper and lower Hubbard bands through the strong correlation effect, the low energy physics is solely described by the coherent 
band near Fermi level and the Hubbard bands appear only as high energy charge fluctuations such as the incoherent absorption 
bands at high frequencies of order 0{U), for example, in optical conductivity. This picture is correct with respect to the charge 
degree of freedom because the DMFT is designed to capture the strong correlation of charge dynamics by being a single site 
impurity model. However, it is physically rather obvious that the localized Hubbard bands can still contribute to low energy 
physics through spin fluctuations. In order to capture this low energy spin degree of freedom, it is, however, necessary to study 
the lattice model - not a small cluster but thermodynamically large lattice. Then we have to give up all the merits of the DMFT. 
At the moment, there is no satisfactory microscopic theory for the lattice model which faithfully treats the strong correlation of 
lai-ge U. 

To this end, we note that there exist two spin correlation functions in our two component model: Xs the one from the local 
spins S(r) and Xf the one from the itinerant spin density s(r) — Ca(r)aorpcp(r). Counting the coupling term to one loop order 
(equivalent to the RPA), the dressed spin correlation functions of the Hamiltonian (1) are written as follows. 



XsHq,^) = Xoi(q.^^)-^'-xoj(q,£i) O) 

X7'(q,£i) = Xo}(q.^^)-/-Xo.5(q,£i) (4) 

where Xo.s is introduced in Eq. (|2|i and Xo.f is the noninteracting spin susceptibility of the conduction band of the fermions. The 
diagrammatic illustration of the derivation of Eqs.(3) and (4) is shown in Fig.l. The noninteracting spin susceptibility Xo.f is 
written as 
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FIG. 1: (a) Interaction vertex of and bare spin correlation functions Xo,s ™d Xo,/- (b) The grapfiic illustration of summations of the 
infinite series for the dressed spin correlation functions and X/ as defined in Eqs.(3) and (4) 
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(5) 



where ^(k) = ■y/e^(k) + A2(k) and the itinerant fermion dispersion e(k) is given by a tight binding model 



e(k) = — 2f(cos(A:t) +cos(A:y)) —2t cosikx) •cos(A;,,) — ^. 



(6) 



For calculations in this paper, we chose t = — 0.4f, and /j = —0.8 k. The overall energy scale t and the choice of parameters 
t , /J will be discussed later with the numerical results. For the superconducting state (SS), we assume a canonical d-wave gap 
function A(k) = Ao[cos(A:r) — cos(^,,)] and for the normal state (NS), we set Aq = in Eq.(5). 

Having two degrees of freedom in the model, two spin susceptibilities and X/ should be calculated on equal footing. 
Previous studies of the local spin correlation embedded in the fermion bath^^ considered only the imaginary part of Xo./ (so 
called Landau damping) in Eq. ^ to damp the spin wave excitations of Eq. dU and the real part of Xo.f is assumed either already 
included in the definition of the bare local spin dynamics described in Eq. or having negligible effects. In fact, when the 
coupling g is weak, this approach is reasonable. But in the strong coupling limit when the dimensionless coupling constant 
X = ■ XQ.f{Q,0) ■ Xo.s(Q,0) ~ <3(1), it is crucial to include both the real and imaginary parts as in the above equations (O 
and dUi. As we can see in the next section, in the strong coupling limit both dressed spin susceptibilities Xs(q,ii) and Xf{q,Q.) 
become a mixture of the local spins and the itinerant fermions and they assimilate to each other with increasing the coupling 
strength X. 

In passing, the Eqs. (3) and (4) are loop expansions (one loop order) but not a coupling constant expansion. Therefore, the 
strong couphng limit of A- ~ (9(1 ) is not a problem, but the higher loop diagrams - for example, vertex corrections in a standard 
many body terminology - need to be worried. This question is the beyond the scope of the current work. 
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FIG. 2: (Color online) (a-b, in the left column) The dressed itinerant spin susceptibility ImXf{({ = {h, 1 /2),Q.) and the bare spin susceptibility 
^'"Xo.y respectively, in the superconducting state. Parameters are A^c = 1-lf, Aq = 0.2t, and X = 0.8; (c-d, in the right column) 

ImXf{q,Q.) and /mXo,/(q,fi), respectively, in the normal state (Aq = 0). 



III. NEUTRON SCATTERING 



In order to study the INS experiments, we calculated the fully dressed dynamic spin susceptibilities Xs(Qj ^) ™d X/(Qi ^) of 
Eqs. (3) and (4). As mentioned in the previous section, in the strong coupling limit of A, ^ 0{1), the behaviors of X5(q,fi) and 
Xf{q,Q.) become qualitatively similar each other Therefore, we conveniently discuss the numerical results of Xf{q,Q.) in this 
paper. But for completeness, we also show the numerical results of Xs(q, ^) as well as the total spin susceptibility Xror —Xf + Xs, 
too. 



A. x/(q,i^) along (0,1/2) ^(1,1/2) 

Figure |2a) shows /mX/(q,n) scanned along q — (/z, 1/2) in the SS. The superconducting gap Aq — 0.2t, the bare spin gap 
AsG = 1.1? (it is not the physical spin gap), and the dimensionless coupling constant A, = 0.8 {g^ = 0.95 eV^) were chosen. 
The main effect of the coupling is to renormalize down the bare spin gap energy A^g below the particle-hole excitation gap of 
Xo./(q,ii) 2Ao), which then forms a sharp resonance peak at Q = (1/2, 1 /2). Centering from this resonance mode, both 
the downward dispersion branch and the upward dispersion branch span out. The origin of the upward dispersion is apparently 
from the local spin wave mode (see Eq. (|2]i) and the origin of the downward dispersion is the itinerant spin excitations of Xo.f- 
The latter fact can be identified in Fig.|2|b) which shows the non-interacting fermion spin susceptibility /mXo,/(q,ii) scanned 
along q — {h, 1/2) in the SS. The shape and strength of the downward whisker like excitations in ImXoj is sensitive to the Fermi 
surface (FS) curvature (controlled by fj and t ), and the size of the d-wave gap A(k) = Ao[cos(^v) — cos(^y)]- 

With the coupling strength X = 0.8, the dressed fermion spin susceptibility X/(q,ii) obtains features of both the local spin 
susceptibility Xo.s and the itinerant spin susceptibility Xo.f- In particular, the high energy parallel branches in Xo.f (see Fig.|2|b) 
and (d)) are overwhelmed by the spin wave like excitations of Xo.s in the dressed susceptibility X/ as seen in Fig. |2a) and 
(c). With a smaller coupling strength (X < 0.5) the two spin susceptibilities Xs(q,ii) and X/(q,ii) retain more of their bare 
characteristics of the spin wave excitations and the itinerant fermion susceptibility, respectively, and the resonance peak at 
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SS NS 



FIG. 3: (Color online) (a-b, in the left column) The dressed local spin susceptibilities /raXs(l = (h, 1 /2),Q.) and the total spin susceptibilities 
ImXiot{(l = {h, 1/2), f2) = ImXs + IiiiXf' respectively, in the superconducting state, (c-d, in the right column) ImXsil,^) and ImXtot{(l,^), 
respectively, in the normal state. Parameters are the same as in Fig.2. 



Q = (1 /2, 1 /2) is not formed. 

Figure|2c,d) are the same plots as in Fig.|2ja,b) but in the NS. First, the resonance peak becomes severely overdamped having 
only a hump like structure in /m%/(q,i2). Second, the low energy downward whisker like dispersion disappears because the 
free fermion susceptibility 5(o,/(q,ii) in the NS (see Fig.|2d)) has no such structure. Lastly, the high energy upward dispersion 
remains almost similar to the case of the SS. The results of Fig.|2a,c) successfully reproduce the main features of recent neutron 
scattering experiments in HTC^~", ie., the resonance mode in SS, the hourglass shape of the upward and downward dispersions, 
and their drastic change between superconducting and normal states. In particular, these results strikingly resemble the INS 
data of YBCOe.e"- However, we need a reservation for applying our result to Li.875B;f=o.i25C04*^, which has a static stripe 
order and extremely low = 2.5 K^^. In particular, the presence of the stripe ordering is likely to change the spin dynamics 
significantly'^"'^ and introduce the a-b plane anisotropy. 

For comparison, we also show the numerical results of the local spin susceptibility 5(5(q,i2) as well as the total spin suscep- 
tibility Xrof(q,ii) = +Xs(Q:^) in Fig- 13 As mentioned above, the overall behavior of 5(s(q,ii) is indistinguishably 
similar to X/(q,ii) and so is Xtot{(l;^)- This is the typical feature of the strong coupling limit of A- ^ 0(1) and remember that the 
value of A, = 0.8 used in our calculations was not an arbitrary choice but was determined by the {k,k) resonance condition. As 
expected, however, a fine difference exists so that, in general, Xs(q,i2), in comparison to Xf{q,Q.), has a slightly more feature of 
the local spin dynamics at higher frequencies and a slightly less feature of the itinerant fermion spin dynamics at low frequencies 
and vice versa. For example, we can see a bit weaker downward dispersion branch in X5(q,£2) in the SS (Fig. 3(a)) than in 
X/(q,n) intheSS (Fig.2(a)). 

Also although we plotted the total spin susceptibility as Xiotai = X/ + Xs, there is an ambiguity about whether the contributions 
from the local spin fluctuations Xs and the itinerant spin fluctuations X/ to the INS measurement should be equal as we tentatively 
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n (in unit of t ) 

FIG. 4: (Color online) (a) Plots of bare susceptibilities i?eX(7s(Q'^)' S^^^Xo.fiQj^)' ™d ^^^™Xo./(Qif^)! respectively in superconducting 
state. Parameters are A^g = l.lf, Aq = 0.2/, and X = 0.8. The vertical dashed line is a guide to the eyes indicating the position of pole at 
Q. = 0.28f. (b) The same as (a) in normal state (ie. Ao = 0). 

assumed here because the form factors of the local and itinerant spins, in principle, should be different. However, our ignorance 
of the relative strength of the %f and fluctuations does not affect our phenomenology because the physical coupling strength 
between two spin fluctuations is determined by the effective dimensionless coupling A,(^) = •Xo./Cq^O) -Xo.sCq^O) and the 
value of A,(Q) is determined once and for all by the (7i,7t) resonance condition and all other physical quantities are calculated 
without further ambiguity. Keeping this point in mind, the comparison of our numerical calculations to the INS experiments 
should be qualitatively the same whether we use the results of X/ or X5, or %tot- 

B. Origin of the resonance mode 

The mechanism of forming the resonance mode in our model is illustrated in Fig. |4] When the inverse of the dressed suscep- 
tibilities of Eq. (O and Eq. (|4|i crosses zero (which occurs simultaneously in both susceptibilities), the dressed susceptibilities 
develop a resonance mode: a bound state or an overdamped mode depending on the presence and strength of the imaginary 
part at the position of pole. In Fig.|4]we plot separately /?eXo ^(Q,i2), ReXQj{Q,D.), and ImXQj{Q,D.) to make this point clear 
Figure HJa) is the case of a SS, where the pole of = Qi^) occurs at Q.,-gs ^ 0.28f. At this frequency the damping from 

ImXo.f is very weak below the p-h excitation gap, so that the pole becomes a sharp resonance peak. FigureHJb) shows the case 
of the NS (Ao = 0) with the same parameters as in Fig.|2c). The position of pole occurs at a little higher frequency (Q. ~ 0.5?) 
compared to the case of the SC phase (Fig.|3a)). But this pole is strongly damped by /m^o./ (green line) that is linearly increas- 
ing with energy, and this linearly increasing damping shifts down the actual position of the peak to Clres ~ 0.35f (the maximum 
height position in Fig.2(c)). We note that this overdamped resonance peak at NS is consistent with the data of Refiii as shown in 
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Fig.Ec). 

The resonance mode found in our model has physically different content than the resonance mode in the Fermi liquid type 
theoriesi^— . The line of /?e%Q^(Q,n) in Fig. |4]is not a simple inverse of a static potential (for example, jjj^ in a RPA 

calculation of Hubbard model as ir>ii"— ) but it carries its own dynamics and spectral density. Therefore, the resonance mode 
formed by coupling of two dynamic susceptibilities Xq j and Xo.s should carry the spectral densities from both the local spin wave 
and the fermion particle-hole continuum. In our coupled two component spin-fermion model, the upward excitation branch and 
the resonance mode appear from a pole of Eq. ([3]) and Eq. (|4| for a given q but the downward excitation branch is made of 
particle-hole excitations of fermions in the d-wave SS and does not constitute a pole in Eq. ([Sj and Eq. This is in contrast 
with the Fermi liquid type theories '^^'^ where both the downward branch and the resonance mode are constructed by the pole 
of a RPA type spin susceptibility. 

C. Constant energy scans 

In the left column of Fig. |5](a-d), we show the constant energy scans of Xf{q,Q.) in the SS for £1 — 0.2f ,0.28f,0.6/, and 
0.8f , respectively. Constant energy scans of neutron scattering data of YBCO^ and LBCO^ show peculiar patterns of IC peak 
positions in {qx,qy) momentum space at different energy cuts. In particular, the 45 deg rotation of the patterns from a low energy 
scan (below the resonance energy Clres) to a high energy scan drew special attention and several theoretical explanations have 
been proposed Results of Fig.|5]demonstrate that the two component spin fermion model can consistently explain this 
phenomena, too. 

Figure|5jb) in the left column is the scan of Xf{q,Q.) at the resonance energy, Q.,es — 0.28f with the same parameters as in 
Fig.|2ja). It shows a very intense peak at (1/2,1/2) indicating a very sharp resonance not only in energy but also in momentum 
space. Fig. |5|c) and Fig. |5jd) are the scans at higher energies than the resonance energy and Fig. |5ja) is a scan of lower energy 
cut. We colored the highest intensity positions with black color to emphasize the clear patterns. The lower energy scan (Fig.|3a)) 
shows the IC peaks at (1/2 ± 5, 1 /2) and (1 /2, 1/2 ± 5) forming a diamond shape pattern. The higher energy scans (Fig. |5jc) 
and Fig.lSfd)) show that the IC peak positions at (1/2 ± 5, 1/2 ± 5) and (1/2 ± 5, 1/2 =F 5) forming a square shape pattern which 
has the symmetry of the 45 deg rotated from the low energy pattern. The results of Fig. |5] excellently reproduce the observed 
patterns of the constant energy scan data of neutron experiments reported in YBCO^ and LBCO^. This is rather surprising for 
LBCO since this compound is known to develop a static stripe ordering and our model has no ingredient for the stripes as we 
mentioned before. 

In our model we can trace the origins of the IC peak patterns. The low energy IC peaks and diamond shape pattern is basically 
a reflection of the band structure and d-wave superconducting gap. The high energy IC peaks and the square shape pattern 
has more complicated origin. At and above the resonance energy the dressed spin susceptibility X/ is the result of a strong 
interplay between the local spin correlation and the itinerant spin correlation. Therefore the high energy scan pattern is the 
result of a subtle interplay/competition between Xq.s{(1,^) and Xo./(q,ii)- The presence of IC peaks at high energies itself is 
the manifestation of the high energy spin wave dispersion spanning from the AFM wave vector Q; so the incommensurability 
increases with energy. However, whether the pattern will be a square or a diamond shape has no universal mechanism. We tested 
various combinations of parameters f , fj, Asg, Aq and X. The low energy diamond shape pattern is robust within our model. As 
to the patterns of higher energy scans, although the square shape is the dominant one, it is not absolutely robust; with different 
parameters the diamond pattern can appears, too. Therefore, we think that the 45 deg rotation of the IC peak patterns may not be 
an universal feature of HTC; it can change with doping and for different cuprate compounds. This non-universality is also seen 
in the scan of Xs{(l,£l) at Clres = 0.6f in the center column of Fig.|5jc). However, as repeatedly emphasized, the patterns of the 
constant energy scans for %5(q,£2), %,o,(q,£2), and X/(q,£i) are basically the same each other in the strong coupling limit. 

D. Parameters of our phenomenology 

To make a comparison of our calculations with experiments, it is important to fix the energy scale of the model. The tight 
binding band of Eq. (|6]l is widely studied to fit the ARPES data and the estimate of t varies from 150 meV to 400 meV depending 
on the doping and different cuprate compounds^^>. Our calculation results are in good agreement with neutron experiments in 
terms of energy scale if we choose t ~ 150-180 meV. This value of t corresponds to the low end of the estimates from ARPES 
experiments. One possible reason for it is that the extraction of t value from ARPES is carried by fitting the whole Brillouin 
zone (BZ) of the quasiparticle (q.p.) dispersions. As a result the high energy dispersion sets the overall energy scale t. However, 
the low energy spin susceptibility is determined by the low energy particle-hole excitations near Fermi level and irrelevant with 
the high energy q.p. excitations. With this reasoning it is quite possible that the effective t value near FS is in fact much reduced 
by a renormalization due to the strong correlation effect. 
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FIG. 5: (Color online) Constant energy scans of /mX/(q,fl) (left column), /mXs(q,fJ) (center column), and /mX/o/(q,fi) (left column) at (a) 
£1 = 0.2t, (b) = 0.28f, (c) £1 = 0.6?, and (d) fl = 0.8?, respectively. In all cases, parameters are Asc = 1. U, Aq = 0.2t, and A. = 0.8. 
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(a) (b) 



FIG. 6: (a) A typical diagram of the Free energy of the model Hamiltonian H of Eq.(l) in one loop approximation, (b) A typical diagram of 
the fermion self energy in one loop approximation calculated in Eq.(7). The diagram (b) can be obtained by cutting any fermion line of the 
diagram (a). 

The degree of incommensurability and the strength of the downward whisker-like dispersion in the fermion susceptibility Xo,/ 
(see Fig.|2b)) is controlled by the FS curvature - which is tuned by t and the chemical potential /j - and the SC gap size Aq; this 
property is true even for the Fermi liquid theories or one component spin-fermion theoriesi^— . These parameters t , /j, and Aq 
should be independently determined by other experiments such as ARPES^^, tunneling^^, etc. Therefore, after fixing the overall 
energy scale of the model by t, the genuinely free fitting parameters of our phenomenological model are only two: the coupling 
strength g (or equivalently X) and the bare spin gap Asg- For all calculations in this paper we used t = — 0.4f, /j = — 0.81f, 
Ao = 0.2t, Asg = l-lf, and A.=0.8. 



IV. ARPES AND HIGH ENERGY KINK 



As a consistent check for our phenomenological model, we calculate the renormalized band dispersion with the same param- 
eter set that we have used for the neutron scattering in previous section. In this calculation, the most important parameter is the 
coupling strength between fermions and spin fluctuations, for which there is no direct experimental measurement nor a reliable 
theoretical estimate from a microscopic Hamiltonian. In our phenomenology, this value X = ■ Xo,/(Q,0) • Xo.s{Q,0) = 0.8 
was determined by the condition to produce {k,k) resonance mode both in the superconducting and normal states. So we can 
crosscheck the consistency of our phenomenology by comparing the outcomes of the q.p. renormalization from the spin-fermion 
interaction to the ARPES experiments. 

The selfenergy of the fermion q.p. is calculated in Born approximation with the fully dressed local spin fluctuations X5(^,(b) 

as 



^ J 71 CO + CO — Zk+q + 'i 

where n(co) and /(co) are the Boson and Fermion distribution functions, respectively. Notice that this calculation of the fermion 
selfenergy is the same one loop approximation as the calculations of the spin susceptibility renormalization in Eqs.(3) and 
(4), guaranteeing the consistency of our phenomenology. This is graphically demonstrated in Fig. |6] We then calculate the 
renormalized q.p. spectral density as A(^,co) = /mG« (^,Co) with G«(^,Co) = m e. m) ' 

In Fig. |7] we show the contour plots showing the intensity of the q.p. spectral densities with and without the selfenergy 
correction. Fig.|2la-b) show the dispersions along the nodal direction (0,0) — > (7i,7i) and Fig.|7|c-d) show the dispersions along 
the near antinodal direction (0,0. 57i) — > (7i,0.57i). Two features are distinctively seen as results of the spin-fermion interaction: 
(1) the overall q.p. dispersion is renormalized by a factor of ^ (9(2). This is consistent with the input X — 0.8 because the wave 
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(h.O.S) in unit of n (h,0.5) in unit of n 



FIG. 7: (Color online) Spectral densities of quasiparticle dispersions: (a, c) Renormalized band dispersions along (h,h) and (h, 0.5) cuts, 
respectively, (b, d) Bare band dispersions along (h,h) and (h, 0.5) cuts, respectively. 

function renormalization factor is Z « (1 + A-); (2) much interesting point is that the q.p.s with a certain energy below the Fermi 
level are so strongly renormalized that the bare dispersion is not continuously renormalized but is detached from the low energy 
dispersion and forms a broad spectral puddle at further high energy region (see Fig.|7ja) and (c)). 

This second feature is commonly observed both for the nodal and antinodal directions regardless of some differences of the 
fine details such as the breaking points of dispersion in energy and momentum space and the intensity of the broad spectral puddle 
at high energies. This behavior is remarkably similar to the so-called high energy kinks ('^ 34QmeV) observed in the ARPES 
experiments with BSCCO and LBCO by T. Valla et af The breaking point of the dispersion occurs around co ~ 1 .5 — 2f in our 
model calculations. In fact, if we assume t ^ 150 — ISOmeV as discussed before, the kink energy we calculated corresponds to 
^ 225 — 3>6QmeV, consistent with the experimental data. Our results even reproduce the overall differences of the dispersion and 
kink behavior between the nodal direction and the antinodal direction as observed in experiments.^^ 

We can trace the origin of this jump or breaking of the q.p. dispersion and found that it is caused by the upper bound of the local 
spin wave excitations X5(q,ii)- The local spin wave excitations, defined in Eqs.(2) and (3), disperses from the lowest energy Ajg 
(renormalized one) at q = Q = (7t,7t) to the highest energy at the magnetic zone corners q = (27t,27t) and its equivalent points. 
Even after dressed by fermions as in Eq.(3) this damped spin wave excitations has a similar upper bound. With the parameters 
of our model, this high energy upper bound of the local spin fluctuations is limited at around 2t. Physical meaning of it is that 
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the local spin excitations exist only up to It and the fermion q.p. cannot be scattered beyond this energy scale. Consequently 
the real part of selfenergy calculated with Eq.(|7]) develops a rapid variation at around 2f in frequencies and q.p. pole is not 
formed beyond this energy scale in the dressed fermion Green's function Gpf^, Co) = ^—^ — -. Hence we can understand 

the origin of the high energy kink from the high energy upper bound of the local spin excitations^. Consistency between our 
model calculations and experimental observation of the high energy kink strongly supports the existence and strength of local 
spin excitations which have a upper energy scale around 35QmeV. In contrast, the itinerant spin fluctuations has a long tail of 
the particle-hole continuum excitations up to the band width 8f) which would be a couple of eV at least. Therefore, the 
one component spin-fermion model with only the itinerant fermions would have a difficulty to explain the high energy kink 
phenomena. 

The reproduction of the high energy kink feature in the renormalized q.p. dispersion also implies that not only the low energy 
spin excitations near (ti, 7i) plays an important role but also the high energy spin excitations can play an important role. However, 
it doesn't necessarily mean that the coupling of the high energy spin excitations to fermions has a comparable strength as the 
coupling of the low energy spin excitations to fermions. In fact, the effective coupling between fermions and the local spins A-(q) 
becomes much weaker with q away from Q = (7t,7t). However, the phase space of the local spin excitations rapidly increases 
with increasing q from Q = (7t,7t), which compensates for the weakness of coupling. 

Within our resolution, we didn't find noticeable features in the low energy q.p. dispersion, i.e., the low energy kink, which 
might be related with the {%,%) resonance mode. Possibly it is because the spectral weight of the resonance mode in NS (see 
Fig. [3jc)) is not sufficiently dominant over the total spectral density of the spin fluctuations spread over the whole momentum 
space, and/or more possibly because the low energy kink is in fact not an abrupt kink but rather a gentle variation of the dispersion 
slope as seen in the recent experiment^. 



V. CONCLUSION 



In summary, we proposed a phenomenological two-component spin fermion model motivated by the neutron scattering exper- 
iments in HTC. With the two spin degrees of freedom of the local spins and the itinerant spins, our calculations of the dynamic 
spin susceptibilities coherently reproduced the essential features of the neutron experiments in HTC: the hourglass dispersions, 
resonance mode, their changes in normal and superconducting states, and the IC peak patterns of constant energy scans. Al- 
though our approach is a phenomenology, considering that there are genuinely only two free fitting parameters, i.e., the coupling 
constant A, and the bare spin gap A^g, the successful reproduction of the several key features of neutron experiments with one 
set of parameters is quite encouraging. Then with the same model parameters, we calculated the renormalized fermion q.p. dis- 
persion and reproduced both the nodal and antinodal high energy kinks with the correct energy scale in agreement with ARPES 
experiments^^. It further strengthened the justification for our phenomenology. 

Finally, the main message of this work with the two component spin-fermion phenomenology is to demonstrate that there are 
compelling experimental evidences^^ for the presence and its important role of the local spin degrees of freedom in addition to the 
fermionic quasiparticles in the cuprates. Interestingly, there are also accumulating experimental evidences for the coexistence of 
the itinerant electrons and local moment of spins in the recently found iron-based superconducting compounds^^'^^"*. The pressing 
question is now what the microscopic theory is for the phenomenological two-component spin fermion model; in other words 
how the local spin degrees of freedom survives doping from the parent insulating cuprate compounds, or more generally how 
the local moments and the itinerant fermions coexist in the strongly correlated metallic systems such as cuprate and pnictide 
compounds. 

The author (Y.B.) was supported by the Grant No. NRF-2010-0009523 and NRF-20 11-0017079 funded by the National 
Research Foundation of Korea. 
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